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1 Introduction 

The process of vector boson production in association with jets is a very important Standard 
Model process at hadron colliders. VF-boson production in association with jets is a typical 
background for many new-Physics searches, and it is also interesting as a process in itself 
to study jets at hadron colliders. The production of a Z-boson in association with jets 
present similar features, since for example when the Z decays in a neutrino pair one can 
have a "missing energy + jets" signature [1], and it is also a reducible background for ZH 
production, similarly to the W + 2 jets case for WH production. Moreover, the production 
of a Z-boson in association with 2 jets is also one of the backgrounds to Higgs-boson 
production via vector-boson fusion (VBF), for example when the Higgs decays in a t~^t~ 
pair. In particular, the signal-background separation between Hjj via VBF and Zjj relies 
on the very different kinematical properties of the jets produced in association with the 
heavy boson [2, 3]. 

For the above reasons, it is important to reach a high-level of accuracy in the predictions 
for the production oi Z + 2 jets at hadron colliders. The full Next-to-Leading Order (NLO) 
corrections for this process have been computed in refs. [4, 5] for the QCD production, 
and in ref. [6] for the electroweak production. Nowadays, one of the methods to improve 
upon existing NLO computations is to merge them with parton showers. This has been 
realized with the MCSNLO [7, 8] and the POWHEG [9, 10] methods for a number of processes 
of increasing complexity [11-16]. In particular, in the recent past, W + 2 and -|- 3 jets 
have been implemented with the MCSNLO method [17, 18], and H -^2 jets using POWHEG [19]. 
Given the role played by jet activity in many searches where Z + 2 jets is a background, it 
is important to achieve the same level of accuracy for this process as well. 

In this paper NLO corrections for the QCD production of Z/^ + 2 jets merged with 
parton showers according to the POWHEG method are presented for the first time. Effects 
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due to photon exchange and ah the spin correlations of the decay products of the Z/7 
intermediate state have been fully taken into account in the implementation, although 
in the following we will refer to the process simply as "Z + 2 jets production". We will 
concentrate here mainly on the implementation procedure and on the discussion of the 
results we obtained, in presence of different cuts, leaving a broader phenomenological study 
for future work, although a comparison with the first available LHC data is presented. 

2 Method and details of the implementation 

To simulate with NLO matched with parton shower (NLO+PS) accuracy the Z + 2 jets 
process, the POWHEG BOX package [20] has been used. The POWHEG BOX is a program that 
automates all the steps described in ref. [10], turning a NLO calculation into a POWHEG 
simulation. The details of how the program works have been largely described in ref. [20] , 
and therefore will not be repeated in this paper. In this section, we first summarize how the 
inputs needed by the package to work were obtained, and we then describe some technical 
details relevant for this implementation. 

The Born matrix elements have been computed using the helicity-amplitude technique 
of refs. [21, 22]. By keeping the amplitude uncontracted with respect to the gluon polariza- 
tion vector, the spin-correlated Born matrix elements {B^^, defined in eq. (2.8) of ref. [10]) 
are easily obtained. The color-connected squared amplitudes {Bij, defined in eq. (2.97) of 
ref. [10]) for this process are not proportional to the Born squared amplitudes. They have 
been computed inserting the color factors correspondent to the product of two color opera- 
tors T.j-Tj while computing the terms needed to build the Born squared matrix elements.^ 
Real corrections have been obtained using MadGraph4 [23], whereas the virtual corrections, 
first computed in ref. [24], have been obtained linking the POWHEG BOX with BlackHat [25], 
using the interface proposed in ref. [26]. A sample of Feynman diagrams that enter the 
Born and the virtual contributions is reported in fig. 1, whereas some numerical values for 
the finite part of the virtual corrections are reported in Appendix A. 

Several checks have been performed at this stage: spin-correlated and color-linked 
squared amplitudes have been checked numerically, by comparing the soft and collinear 
limits of the real matrix elements with the expected value obtained from the factoriza- 
tion formulas for collinear and soft emissions, respectively. Thanks to the aforementioned 
interface, BlackHat also returns the coefficients of the double and the single pole of the 
one-loop corrections: this allows for an extra cross check of the color-linked amplitudes, 
or, viceversa, a check that the interface for the 2 codes is working properly. As a final 
important check, several distributions have also been compared against NLO predictions 
obtained from n-tuples generated with BlackHat+SHERPA [27, 28], and agreement has been 
found. 

^Similarly, the relative weights of each colour structure present in the Born processes, in the limit 
of large number of colours, have been computed combining the terms needed to build the Born matrix 
elements. These weights are used to probabilistically assign a planar colour structure to the underlying- 
Born kinematics, from which a color flow is attached to the generated events, using the prescription described 
in refs. [10, 20]. 
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Figure 1. Sample graphs for the Born and the virtual contributions to the Z/7 + 2 jets production 
process. 

The significant property of a process like V + 2 jets, that makes it more difficult to 
implement than 1/+ 1 jet, is the fact that the Born matrix elements are singular in several 
regions. In order to suppress the singularities of the underlying-Born amplitudes, instead 
of using sharp generation cuts, we used a suppression factor that vanishes when at least one 
of these regions is approached. As explained in [11], this means that the underlying-Born 
kinematics is generated according to a modified B function 



where B is the inclusive NLO cross section at fixed underlying-Born variables, and, for the 
case at hand, we used 



light partons in the underlying-Born kinematics and si^2 is the invariant mass squared 
obtained from their momenta. Since the Born-like kinematics, from which the radiation 
is generated, is distributed according to eq. (2.1), the generated events will be assigned a 
weight proportional to 1/F. To obtain the results shown in the next section, we have set 
kis = kps = 2, Aprj, = 10 GeV and Am = 5 GeV. We also notice here that a similar method 
has been recently used in ref. [19] to implement in POWHEG the H + 2 jets process. 



(2.1) 




(2.2) 



is the Born phase space, pT,i and pT,2 are the transverse momenta of the 2 outgoing 
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3 Results 



3.1 Comparison with NLO predictions 

In this section we present results obtained after showering and hadronizing the partonic 
events generated with POWHEG. We have used PYTHIA 6.4.25 with the "Perugia 0" tune 
enabled. We have considered Z(-^ e^e~) + 2 jets production at the LHC, with an hadronic 
center-of-mass energy ^/S = 7 TeV. 

The values of physical parameters entering the computation are the following: 

mz = 91.1876 GeV, rz = 2.49GeV, a"^ = 128.802 , sin^ 6*^/ = 0.23 , (3.1) 

and we have used the CTEQ6M [29] parton distribution functions. In the computation 
of the B function, the renormalization and factorization scales have been chosen equal to 
Ht/'^, where 

Ht = \Jm%+ pip z + PT,i + PT,2 , (3.2) 

Pt,z is the transverse momentum of the Z-boson, and all the quantities are computed using 
the underlying-Born kinematics. 

As far as technical parameters or special options of the POWHEG BOX program are con- 
cerned, we have kept them equal to their default value, and we have not used any folding 
in the integration of the radiation variables [11, 20]. In so doing, the fraction of negative- 
weight events amounts to be 21%. These events were kept in the final sample used for the 
analysis, although of course not all of them will contribute to the plots, since some will not 
pass the cuts. 

Jets have been defined according to the anti-fcx algorithm [30-32], setting R = 0.4, 
and the following cuts have been enforced in the analysis: 

66.328 GeV < m^+^^ < 116.048 GeV , pT,e > 20 GeV , jyej < 2.5 , 

PTj> 30 GeV, \rij\<4A. (3.3) 

Events are accepted only if there are at least 2 jets passing the above cuts. We also remind 
that we have not used isolation cuts for the leptons, and that photon radiation off leptons 
have been switched off in the shower,^ in order to allow for a clear comparison with the 
fixed order result. 

In fig. 2 we show the transverse momentum spectra of the reconstructed Z boson 
{pt,z), of the positron {pT,e+)^ ^-iid of the hardest and second-hardest jet (pt, jn Pt, j2 ) ' 
whereas in fig. 3 the invariant mass of the system made by the 2 hardest jets (mjj) and 
their azimuthal separation (A^jj) is shown. 

The agreement between the NLO predictions and the results after the shower and the 
hadronization stage is very good for these inclusive observables. Shapes are essentially 
unmodified, and very small diff'erences can be observed in the lower panel of each plot, 
where the ratio between the showered result and the NLO is plotted.^ We have checked 

^Photon radiation has been switched off by setting MSTJ(41)=3 in PYTHIA. 

^Errors in the lower panels of all the plots in this article represent the error on the ratio between two 
results, obtained propagating the errors of each set. 
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that they are of the same order of the theoretical uncertainty of the NLO result, obtained 
by changing the renormalization and factorization scales by the usual factors {1/2,2}. It 
is also worth noticing that effects of the same order have been observed in other NLO+PS 
implementations of similar complexity [17, 18], although in presence of different cuts. 
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Figure 2. Comparisons between POWHEG (interfaced to PYTHIA) and the NLO results at the LHC 
pp collider {^/S ~ 7 TeV), for the transverse momenta of the reconstructed Z-boson, of the positron 
(upper panel), of the hardest and of the second- hardest jet (lower panel). Vertical bars correspond 
to statistical errors. 

In the upper panel of fig. 4 we show instead the transverse momentum of the third 
hardest jet and Htj, which is defined to be the scalar sum of all the jets transverse 
momenta: good agreement is found between the NLO result and the POWHEG+PYTHIA 
prediction for these observables as well. In particular, for values of pr.js far from the low 
Pt region, the results after the showering and hadronization stage agree with the fixed 
order result (which has leading-order accuracy for this observable). In the small transverse 
momentum region, the NLO+PS prediction for ptjs is expected to show Sudakov damping, 
as opposite to the NLO result, which diverges: this is noticeable in the lower panel of fig. 4, 
where we plot PT,j3 removing the lower cut on the transverse momentum of all but the two 
hardest jets. The results obtained just after the parton-shower stage performed by PYTHIA 
are also shown for this observable. Sudakov damping is clearly present in all the curves 
where resummation of soft-collinear emissions is performed, although in the very 1ow-Pt 
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Figure 3. Comparisons between POWHEG (interfaced to PYTHIA) and the NLO results at tlic LHC 
pp collider {\/S ~ 7 TeV), for the invariant mass of the two hardest jets and for their azimuthal 
separation. Vertical bars correspond to statistical errors. 

region non-perturbative effects seem to play a small but noticeable role. 

We have also performed the same analysis using as factorization and renormalization 
scale the Z-boson transverse mass, computed with the underlying-Born kinematics (// = 
Et,z = Y^rfi^ + and we have found similar results. 

As mentioned in the introduction, Zjj production is a background for Higgs-boson 
production via VBF, when the H ^ tt decay channel is considered in the analysis. For 
this reason it is interesting to check how typical distributions look like for the QCD Zjj 
background, in presence of a minimal set of VBF cuts. We chose 

PT,f>20GeV, |2/£|<2.5, 

|r?, I < 5.0 , pT,j > 20 GeV , ptj,^^ > 30 GeV , 

I'^ii -^i2l > 4.0, ??ji-%2<0, 

min (r/jj , r]j^ ) + 0.4 < /f- < max {rjj^ , r/^^ ) - 0.4 , (3.4) 

and we restricted the invariant mass of the lepton pair to lie in the interval [66.328, 116.048] 
GeV. 

The "tagging jets" are the two hardest jets, and we have used as before the anti- 
algorithm, with R = 0.4. For this analysis we have generated the events using as 
factorization and renormalization scale the Z-boson transverse mass, computed with the 
underlying-Born kinematics. Moreover, multiple parton interactions have been switched 
off. 

In the upper panel of fig. 5 we show pt,z and PT,ji, whereas in the lower panel we 
show the azimuthal decorrelation between the two hardest jets and the shifted rapidity 
y^'^^ = Vjs ~ [Vji + yj2)/2) which is an useful quantity to measure the distance between the 
tagging jets and the third hardest jet. 

Since we are using cuts designed to suppress the Zjj background, the statistical sig- 
nificance of the plots in fig. 5 is not optimal. However, we found that the results are not 
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Figure 4. Comparisons between POWHEG (interfaced to PYTHIA) and the NLO results at the LHC 
pp colUder {^/S = 7 TeV), for the tliird hardest jets transverse momentum and for i?rj (upper 
panel). In the lower panel the low cut on the jets transverse momenta is kept only for the hardest 
and second-hardest jet. Vertical bars correspond to statistical errors. 

changed sizeably when going from a NLO to a NLO+PS prediction. In particular y^'^^ is 
peaked at 0, whereas, as a consequence of the VBF cuts, the rapidity of the two tagging 
jets (not shown) is peaked at \yjt^J — 2.5. This is the expected behaviour, since the QCD 
production of Zjj is a process dominated by quark or gluon exchange in the t-channel, 
and therefore one expects to observe a third jet in the central rapidity region, as opposite 
to what happens in processes dominated by a t-channel exchange of a colourless state, like 
Higgs-boson production via VBF or Zjj EW production. In the latter cases, the extra jet 
activity tends to be close to one of the two tagging jets, and, as a consequence, y^'^^ peaks at 
high rapidities and exhibits a dip at y'^'^^ = [33, 34]. The azimuthal decorrelation between 
the tagging jets is also particularly relevant: indeed it was shown that properties of the 
Higgs-boson such as its parity could be measured when it is produced via VBF, by looking 
at the azimuthal separation between the tagging jets [35]. This distribution appears to be 
stable as well for the Zjj background, when going from NLO to NLO+PS accuracy. 
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Figure 5. Comparisons between POWHEG (interfaced to PYTHIA) and the NLO results at the LHC 
pp coUidcr (a/S" — 7 TeV) in presence of the VBF cuts of eq. (3.4): in the upper panel the transverse 
momenta of the reconstructed .Z-boson and of the hardest jet are shown, whereas in the lower panel 
the azimuthal decorrelation between the tagging jets and y"^^ = yj^ — {yj^ + yj^)/2 are reported. 
Vertical bars correspond to statistical errors. 

3.2 Comparison with LHC data 

The ATLAS Collaboration has published a study on the production of a Z-boson in associ- 
ation with jets [36], and good agreement has been found when comparing the experimental 
results with QCD NLO perturbative predictions and predictions from Monte Carlo gener- 
ators implementing LO matrix elements matched with parton showers. In this section we 
show a comparison of our NLO+PS predictions with ATLAS data. The set of cuts we used 
is very similar to that of eq. (3.3). Here we ask for 

66 GeV < m^+g- < 116 GeV , 

PT,e > 20 GeV , |?/e| < 2.5 , Ai?j- e > 0.5 , 

PT,i>30GeV, |?/j|<4.4, (3.5) 

and jets are built with the anti-fcx algorithm, with R = 0.4. 

In fig. 6 we show our predictions, obtained with the sample generated with /z = Ht/2, 
together with the ATLAS data. In the upper panel we show the transverse momentum of 
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the second-hardest jet and the invariant mass of the two hardest jets, for events with at 
least 2 jets. In the lower panel instead, the azimuthal separation A^jj between the two 

hardest jets and their distance ARjj = \J (<I>jj — ^j^)^ + iVji ~ ViT^ the (j) — y plane is 
reported. The outcome of these comparisons is that the agreement between experimental 
data and the POWHEG results is very good. 
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Figure 6. Comparisons between POWHEG (interfaced to PYTHIA) and the ATLAS data {\fS = 7 
TeV), in presence of the cuts of eq. (3.5), for the transverse momentum of the second-hardest jet, the 
invariant mass of the two hardest jets (upper panel), the azimuthal separation and the Ai? distance 
between the two hardest jets (lower panel). Vertical bars on the POWHEG results correspond to 
statistical errors, whereas bars on ATLAS data represent the sum (in quadrature) of the statistical 
and systematical errors reported in ref. [36]. 



4 Conclusions 

In this paper we have described the first implementation of Z + 2 jets production at 
next-to-leading order in QCD, in the POWHEG framework. We have used the POWHEG BOX 
package, which is a program that automates the algorithm first proposed in ref. [9] and 
then described in detail in ref. [10]. The main purpose of this paper was to show that 
no particular problems occur in the implementation of complicated processes with POWHEG. 
This has been recently observed in ref. [19] , and we obtained similar findings. We regulated 
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the divergent underlying-Born process by using a damped B function, and we have shown 
that results are in good agreement with the fixed order for observables where they are 
expected to. Sudakov damping is instead present for observables sensitive to multiple soft- 
collinear emissions, such as the low-px region of the transverse momentum of the third 
hardest jet. 

An important phenomenological application which can be performed having this pro- 
cesses available with NLO+PS accuracy is a study of the QCD radiation patterns of Zjj 
with respect to those present in Hjj production, especially in presence of cuts used in 
VBF searches. We have given an example of such a study by checking that typical ob- 
servables in Zjj production in presence of VBF cuts have the behaviour expected from 
general properties of QCD, and they are not affected sizeably when going from NLO to 
a NLO+PS description. Further studies along these lines could also be performed using 
NLO+PS implementations of Hjj via gluon [19] and vector-boson fusion [34, 37] and tt(+ 
jets) [13, 14, 38], together with the Zjj implementation presented here. 

Moreover, it could also be interesting to merge NLO+PS predictions for Z + 2 jets 
together with the same predictions for lower multiplicities, using for example the approach 
proposed in ref. [39]. These studies are left for future work. 

We have also compared our NLO+PS results with recent ATLAS data, and we have 
found good agreement. 

The computer code for this POWHEG implementation will soon be available within the 
public branch of the POWHEG BOX package. It will be possible to link it with public one-loop 
codes to evaluate the virtual corrections. 
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A One-loop virtual amplitudes 

In this appendix, we report the numerical values of the finite part of virtual corrections for 
the following phase space point: 

= (138.4784456174; 0.0000000000, 0.0000000000, 138.4784456174) GeV, 
Pq = (138.4784456174; 0.0000000000, 0.0000000000, -138.4784456174) GeV, 
Pe- = (39.8695145736; -30.2118421708, 7.3821554172, -24.9464740270) GeV, 

= (61.7949044527; 60.7995170788, -7.3821554172, -8.2178294395) GeV, 
pi = (104.9831601822; -38.2721470356, 44.5354030781, 87.0247353101) GeV, 
P2 = (70.3093120261; 7.6844721276, -44.5354030781, -53.8604318436) GeV. 

(A.l) 

The finite parts of the interference between the tree-level and the one-loop amplitude 
(as returned by BlackHat), computed in the 't Hooft-Veltman scheme and stripped off of 
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subprocess 


/I tree ^ 


2 3ft{(^t'^^^)* -^1-1°°?} 


uu — >■ e~e'^gg 


1.0489864146 E-005 


1.7673954550 E-004 


dd — )• e~e^gg 


1.3620189441 E-005 


2.4684280244 E-004 


uu e~e~^uu 


6.9320177545 E-006 


2.1474734072 E-004 


uu — )• e~e~^cc 


1.8468108910 E-007 


1.8007297663 E-006 


dd — )• e~e~^dd 


1.0023268662 E-005 


3.1647659850 E-004 


dd — )• e~e~^ss 


2.3673784898 E-007 


2.4179197657 E-006 


uu — )• e~e^dd 


6.6505600059 E-007 


1.3323314116 E-006 



Table 1. Partonic subprocesses and corresponding values for the Born squared amplitudes and the 
interference between the tree-level and the one-loop amplitudes. Momenta are given in eq. (A.l), 
whereas physical parameters are reported in the text. Processes that can be obtained by crossing 
are not reported. 



the usual factor 

are reported in Table 1, together with the corresponding partonic subprocesses and the 
Born squared amplitudes. The physical parameters chosen are as reported in section 3.1, 
but here the renormalization scale have been set equal to mz- The results for the virtual 
amplitudes are divided by as/^-K, i.e. they do not contain the extra power of the strong 
coupling constant. We stress that with the BlackHat version that has been used in this 
work these amplitudes are computed in the large nrit limit, and the axial contributions due 
to fermionic loops are not included. It has been explicitly checked that the above results 
agree with the same amplitudes obtained with MCFM, when the aforementioned contributions 
are turned off. Agreement between BlackHat and MCFM has been obtained also using Ht/'^ 
or the Z-boson transverse mass as renormalization scale. 
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